Transient energy excitation in shortcuts to adiabaticity for the time dependent 

harmonic oscillator 
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There is recently a surge of interest to cut down the time it takes to change the state of a quantum 
system adiabatically. We study for the time-dependent harmonic oscillator the transient energy 
excitation in speed-up processes designed to reproduce the initial populations at some predetermined 
final frequency and time, providing lower bounds and examples. Implications for the limits imposed 
to the process times and for the principle of unattainability of the absolute zero, in a single expansion 
or in quantum refrigerator cycles, are drawn. 
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I. INTRODUCTION 

Adiabatic processes in quantum systems are frequently 
useful to drive or prepare states in a robust and con- 
trollable manner, and have also been proposed to solve 
complicated computational problems, but they are, by 
definition, slow. (The definition of "adiabatic process" 
here is the usual one in quantum mechanics, namely, 
a slow change of Hamiltonian parameters keeping the 
populations of the instantaneous eigenstates constant all 
along.) Thus a natural objective is to cut down the 
time to arrive at the same final state, possibly up to 
phase factors, in other words, to find "shortcuts to adi- 
abaticity", by designing optimal adiabatic pathways, or 
by admitting that the populations may not be preserved 
at intermediate times. Several works have recently pro- 
posed different ways to achieve this goal for general or 
specific cases jll— TlOj| . One of the early applications con- 
sidered has beenparticle transport without vibrational 
heating [I|, H, H, 0, Qlj ■ Another important case is fric- 
tionless harmonic trap compressions or expansions for 
state preparation |^-|8|, |ll|, which were first addressed 
with "bang-bang" (piecewise constant frequency) meth- 
ods [ll|; other route is to design by inverse engineering 
techniques a time dependent frequency for which the ex- 
panding modes associated with Lewis-Riesenfeld invari- 
ants [l2[ take the state from the initial to the final po- 
tential configuration without transitions 0, 0] ; this has 
been extended to the Gross-Pitaevskii equation with a 
variational ansatz [J], and has been also implemented 
experimentally to decompress 87 Rb cold atoms in a har- 
monic magnetic trap [13| . In the same vein, Berry has 
provided an algorithm to construct a Hamiltonian H(t) 
for which the adiabatic approximation for the state evolu- 
tion under a time-dependent reference Hamiltonian H(t) 
becomes the exact dynamics with H{t). This algorithm 
has been applied to spins in magnetic fields [H , harmonic 
oscillators [6( , or to speed up adiabatic state-preparation 
methods such as Rapid Adiabatic Passage (RAP), Stim- 
ulated Rapid Adiabatic Passage (STIRAP) and its vari- 
ants Q. Also, Masuda and Nakamura have adapted for 



adiabatic processes @ a ( "fast- forward" ) scaling tech- 
nique to speed up the state dynamics [l4j, with appli- 
cation examples to particle transport or time dependent 
harmonic potentials [5] , and spins or charged particles in 
electro-magnetic fields [T(| . Finally, optimal control the- 
ory has also been used for nonadiabatic cooling under im- 
posed costs [ll|, HH ; and Lyapunov control methods have 
been proposed to speed up quantum adiabatic comput- 
ing without information of the Hamiltonian eigenstates 

In this paper we shall examine the energy "cost" of 
such processes; more precisely, their transient excitation 
energies. Our central study case is the expansion (or 
compression) of a harmonic oscillator, a basic model for 
many operations in any cold atoms laboratory [7|. Intu- 
itively, one expects the transient system energy and the 
time of the process to be "conjugate", i.e., an increase 
of the former when decreasing the later, but the details 
of this relation, and the role played by other parameters 
defining the process (such as initial and final frequencies) 
have to be clarified both for fundamental reasons and for 
the applications. In particular, the energy excitation will 
set limits to the possible speed-up. In a trap which is 
harmonic near the ground state but not for higher ener- 
gies, large transient energies will imply perturbing effects 
of anharmonicities and thus undesired excitations of the 
final state, or even atom loss. The transient excitation 
energy has also implications for quantifying the principle 
of unnatainability of zero temperature, first enunciated 
by Nernst [ljj ■ Fowler and Guggenheim [l8| formulate it 
as follows: "/£ is impossible by any procedure no matter 
how idealized to reduce the temperature of any system to 
the absolute zero in a finite number of operations." They 
identify it with the third law of thermodynamics although 
this is sometimes disp uted. More recently, Kosloff and 
coworkers [II], EH l20j | have restated the unattainability 
principle as the vanishing of the cooling rate in quantum 
refrigerators when the temperature of the cold bath ap- 
proaches zero, and quantify it by the scaling law relating 
cooling rate and cold bath temperature. We shall exam- 
ine the consequences of the transient energy excitation on 
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the unattainability principle at two levels, namely, for a 
single, isolated expansion, and considering the expansion 
as one of the branches of a quantum refrigerator cycle. 

When describing these cycles and indeed in many inter- 
section areas between quantum mechanics and thermo- 
dynamics one finds the need to use the word "adiabatic" 
in two different ways: the thermodynamical one (mean- 
ing that there is no heat transfer between system and 
environment) and the quantum one. Many authors have 
pointed out this duality as an unfortunate source of con- 
fusion. It may prove useful to distinguish them and avoid 
ambiguities and the hassle of detailed explanations with 
a shorthand notation. Following the example of Dirac's 
q- number versus c-number distinction, we propose to re- 
fer to a process as "t-adiabatic" if it is thermodynami- 
cally adiabatic, and as "g-adiabatic" if it is a quantum 
mechanically adiabatic (i.e. slow) process. 



II. BANG-BANG METHODS 

The Hamiltonian for a particle with mass m in a time- 
dependent harmonic oscillator is given by 



H = p 2 /2m + muj 2 (t)q 2 /2. 



(1) 



Let us assume an expansion (compressions are treated 
similarly) with initial angular frequency ujq = uj(0) at 
time t = and final frequency ujf = w(t/) < at 
time tf. 

In the "bang-bang" approach the frequency is shaped 
as a stepwise constant function of time, choosing the step 
values and durations so as to preserve the initial state 
populations in the final configuration. For real trap in- 
termediate frequencies this requires a minimal total ex- 
pansion time 



y/1 - UJf/ujQ 



(2) 



The limit can be realized by only three jumps, i.e., two 
real intermediate frequencies (it cannot be improved by 
using more intermediate frequencies) , but one of the in- 
termediate frequencies should be infinite. 

Up to a constant factor the main dependence in the 
bound © already appears in a simpler process that re- 
produces for (ujf,tf) the initial populations with just 
one intermediate frequency, the geometric average uj\ = 
(wqW/) 1 / 2 , and a total time 



2^/ojQUjf ' 



(3) 



which is a quarter of the corresponding period plj . For 
an initial n-th state of the oscillator the instantaneous 
mean energy during the transient period becomes the 
arithmetic mean of the initial and final energies, 



(H) = h 



LJQ +OJf 



(4) 



Eq. © and the bound © are relevant because the t- 
adiabatic expansion is actually the speed bottleneck in 
quantum Otto refrigerator cycles which use particles in 
a harmonic oscillator as the "working medium" [IJ, [2(| ■ 
The cooling rate R as the cold bath temperature T c ap- 
proaches zero is dominated by the expansion time and 
scales as R oc uif/tf. Since u>f « T c as T c -> 0, the 
dependence of tf on ojf quantifies the unattainability 
principle. In particular, g-adiabatic expansions lead to 
R oc T c 3 scalings, in contrast with the R oc T^ 2 scaling 
achieved with the times in Eqs. © or © (llll20j. 

In Q it was demonstrated, however, that the minimal 
time in Eq. © can be beaten with bang-bang methods 
and inverse engineering methods, see also the next sec- 
tion, by allowing for imaginary intermediate frequencies, 
i.e., transients in which the harmonic oscillator becomes 
a parabolic repeller. It was pointed out [7], that this 
new freedom leads to the absence, at least in principle, 
of a lower bound for the expansion time, which could 
obviously affect the optimal scaling of cooling rates. We 
shall analyze the impact of these ultra-fast expansions on 
the third law in the following sections. 



III. ENERGY BOUNDS FOR 
INVERSE-ENGINEERED TIME-DEPENDENT 
HARMONIC OSCILLATORS 

A. Bound for time-averaged energy 

In this subsection, we will set a lower bound for the 
time-averaged energy in the transitionless expansions and 
compressions of the time-dependent harmonic oscillator. 
A shortcut to adiabaticity taking the n-th state of the 
initial trap to the final n-th state of the final trap up to 
phase factors is achieved [7] by designing the frequency 
from the Ermakov equation 



b + L0 2 {t)b 



6 3 ' 



(5) 



where b is an engineered scaling function which satisfies 
the following boundary conditions at t = and tf, 



6(0) = 1, 6(0) = 0, 6(0) = 0, 
&(«/)= 7>&(*/) = 0,6(i/)=0. 



(6) 



Here 7 = WtJoJuif, and the single and double dots denote 
first and second derivatives with respect to time. The 
simplest choice for interpolating b(t) between and tf 
is a polynomial form, b(t) = 6(7 — l)s 5 — 15(7 ~ l)s 4 + 
10(7 — l)s 3 + 1, where s = t/tf. In this manner the 
n-th stationary state of the initial oscillator will evolve 
according to the "expanding mode" 



/4e _i(n+l/2)/*rft'^ 



(2™n!6) 1 /2 



x e 



T~Ln 



\~JTJ b 



(7) 
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where 1-L n is a Hermite polynomial, and will become even- 
tually, up to a phase, the n-th eigenstate of the final trap 
at tf. At intermediate times |\& n ) does not coincide in 
general with the instantaneous eigenvectors |n) of H(t), 
H(t)\n(t))=e n (t)\n(t)). 

For the n-th expanding mode, the instantaneous aver- 
age energy E n (t) = (V n \H(t)\V n ) is 



E n (t) 



(2n + l)fi 



u: 2 (t)b 2 



b 2 



(8) 



which is in general different, except at initial and final 
times, from e„. The time average of E n is defined by 



1 



E n (t)dt. 



(9) 



To find a lower bound for E n we substitute Eq. (jSJ into 
Eq. ^ , and integrate by parts making use of the bound- 
ary conditions (JB]), 



E n 



(2n + l)fi 
2ujotf 



b 2 



dt. 



(10) 



The integrand has the form of the Lagrangian of a par- 
ticle in an attractive inverse square potential, but the 
minimization problem, i.e., finding an optimal function 
b(t) subjected to the boundary conditions © cannot be 
solved with an ordinary Euler-Lagrange equation since 
there are too many boundary conditions which affect not 
only b but also b and b at the edges of the time interval. 
We can nevertheless find easily, using the Euler-Lagrange 
equation, the quasi-optimal "trajectory" b(t) that mini- 
mizes the integral subjected only to the boundary values 
of b, that is, b(0) = 1 and b(tf) — 7. Since these two con- 
ditions define a broader set of functions than the ones sat- 
isfying ([5]), the quasi-optimal b provides at least a lower 
bound for the time-averaged energy. For the function 



f(t,b,b)=b 2 +u J 2 /b 2 



(11) 



the Euler-Lagrange differential equation /& — 4zfi = is 



b A b = -lo 2 . 

The solution satisfying 6(0) = 1 and b(tf) — 7 is 



b(t) = J(B 2 



2Bs + 1, 



(12) 



(13) 



7 2 + 



uj 2 t 2 and the positive root should 



where B = — 1 

be taken. Substituting Eq. (TTBl into the integral (|TU| . 
we finally obtain a lower bound for the time-averaged 
energy, 

(2n+ l)h 



B n = 



-{(B 2 -oj 2 t 2 )-2u Q t f 



arctanh 



arctanh 



B 2 +B 



1 , 2 f 2 



U tf 



B 

UJ tf 



(14) 



such that E n > B n . When the final frequency uif is small 
enough to satisfy tf <C and 7>1, the lower 

bound has the following simple asymptotic form, 



(2n+l)h 
2u f t 2 f 



(15) 



Incidentally, B n sets also a lower bound for the maximum 
of the instantaneous energy E n (t) . 
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FIG. 1: (Color online) Lower bound for time-averaged energy 
(in units of Eq(0) — fojo/2) as a function oft/ and w/, where 
cjo = 2tt x 250 Hz. 

In Fig. [TJ we plot Bq as a function of t f and oj f for 
ujq = 2tt x 250 Hz, which will be also the initial frequency 
in the following examples. The important point is that 
the transient energy increases not only with decreasing 
final time t / but also with decreasing final frequency to / . 
Figure [2] shows the exponents of the scaling for the bound 
(fT4)) . its asymptotic form (fl"5j) . and the time-averaged en- 
ergy for a polynomial b. In all cases Eq, or Bo oc l/(u>ft 2 ) 
asymptotically. 

A consequence of Eq. dT51) is 



t f > 



(2n + l)h 
2cj f E^ 



(16) 



The interest of Eq. ([TBI) compared to Eq. © is that 
in principle, for fixed ujq, it is possible to beat the bang- 
bang minimal time, but the price is an increase in the 
transient energy. In practice this energy cannot be ar- 
bitrarily large, if only because there are no perfect har- 
monic oscillators. In particular, if we consider in Eq. f| 16[) 
that E n is limited by some maximal value, because of an- 
harmonicities or a finite trap depth, the obtained scaling 
is fundamentally the same as for bang-bang methods, and 
leads to a cooling rate R oc T^ 2 in an inverse quantum 
Otto cycle, although an opportunity is offered to improve 
the proportionality factor by increasing the allowed E n - 
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FIG. 2: (Color online) Dependence of the time-average energy 
for the ground state Eo/Eo(0) on the (a) short time tf ((jJf = 
2tv x 0.25 Hz) and (b) final frequency cjf (tf = 2 ms). In both 
cases Ldo — 2tt x 250 Hz. Bound given by Eq. (|14[) (solid red 
line); asymptotic expression Eq. (|15l) (dot-dashed magenta 
line), and time-average energy for polynomial b (dotted blue 
line). The 8's are the asymptotic exponents of tf and uif 
respectively as they go to zero. 



Independently of the participation of the harmonic 
trap expansion as a branch in a refrigerator cycle, we 
may see rather directly the effect of the above analy- 
sis on a single expansion by assuming that the initial 
and final states are described by canonical density op- 
erators characterized by temperatures To and Tf, re- 
lated by Tf = (uif/uo)To for a population-preserving pro- 
cess. Within the idealized but specific context of a pure 
parabolic potential expansion, the unattainability of a 
zero temperature can be reformulated microscopically as 
follows: The transient excitation energy becomes infinite 
for any population-preserving and finite-time process if 
the final temperature is zero (which requires ujf = 0). 
This excitation energy has to be provided by an exter- 
nal device, so the absence of a lower process time limit 
should not lead us astray, since there remains a funda- 
mental obstruction to reach Tf — in a finite time, in 
the form of the need for a source of infinite power. 



B. Minimization of time-averaged energy 

In order to minimize the time-average energy and ap- 
proach the lower bound, we can use the quasi-optimal 
b{t), Eq. (JT3J), in a central time segment [r, 1 — r], and 
match it at the extremes with two "cap polynomials", 
each of them satisfying three of the boundary conditions 
at Eq. © (at t = or tf), plus three boundary condi- 
tions for 6, b and b at the matching point. The idea is 
illustrated in Fig. [3l 
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FIG. 3: Example of hybrid b combining the quasi-optimal tra- 
jectory in the central segment [r, 1 — r] and "cap polynomials" 
with the right boundary conditions, where ujo = 2n x 250 Hz, 
(jjf — 2ir x 0.25 Hz, tf =2 ms, r = 0.4 (solid red line), quasi- 
optimal square- root b of Eq. ()13[) (dotted black line), and 
polynomial trajectory (blue dashed line). 

The resulting hybrid b takes the form 

(0 < s < t) 



b = 



J2j=Q C 3 S ° 

(B 2 - u 2 t 2 f )s 2 



2Bs + 1 (t<s<1-t) (17) 
(1 -r < s < 1) 



where the coefficients {cj} and {dj} have lengthy expres- 
sions but are easily obtained from the matching condi- 
tions so we omit their explicit forms here. 

Figure 0] demonstrates that this hybrid b can indeed 
minimize the time-averaged energy by making r smaller 
and smaller, approaching the lower energy bound as r — > 
0. A detailed calculation shows that the contribution 
from the caps does not vanish as r — > 0, so the value 
of the time-average bound is reached at the price of a 
singular instantaneous energy, see Fig. 01 



IV. ENERGY VARIANCE AND 
ANANDAN-AHARONOV RELATION 

We shall now discuss the impact of the shortcuts to 
adiabaticity on the standard deviation of the energy 
AH = ((H 2 ) - (H) 2 ) 1 / 2 . This is important because, 
a small averaged energy could in principle be spoiled by 
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FIG. 4: (Color online) Time averaged energy versus r for 
polynomial b (solid blue line), bound in Eq. (|14[1 (red dashed 
line), and hybrid b (dot-dashed magenta line). The contribu- 
tion from the "cap polynomials" (orange line with "*") and 
the central segment (black line with "•") are also depicted. 
Parameters are the same as in Fig. 3. 



a large standard deviation. Anandan and Aharonov [22[ 
found a relation between the time average of the stan- 
dard deviation of the energy and the time of a process 
connecting two given states, irrespective of the Hamil- 
tonian used to connect them. The so called "Anandan- 
Aharonov" (AA) relation provides a lower bound for the 
average uncertainty of the energy, which is extensively 
used to minimize the time tf required for the evolution 
between the two orthogonal quantum states. Based on 
the Fubini-Study metric, the following distance may be 
defined, 



S = 2 



AH(t) 



dt > S Q , 



(18) 



where the minimal value So = 2 arccos(| (^(t = 0)\^(t = 
tf)}\) corresponds to the "geodesic" [2J|. For orthogonal 
states So = tt, and 



where 



AHtf > -, 



_ = £' AH(t)dt 



(19) 



(20) 



but more generally, for arbitrary (possibly non- 
orthogonal) initial and final states, the AA relation is 



AHtf > 



hS 

47T ' 



(21) 



This may be applied to the harmonic oscillator and in 
particular to any process taking an n-th initial eigenstate 
to an n-th final eigenstate, up to phase factors. For the 
ground state we find, see ([7j), 



cos 



'So 



K* (0)1*0 (*/))! 2 



2^/uJoTJJ 



uj - 



(22) 




1.0 1.5 2.0 
logiot^oM 



3.0 



c 



a i 



o u 



8=-l 


(b) 









loglOL^/M)] 



FIG. 5: (Color online) Dependence of the time-averaged un- 
certainty of energy AHo on the (a) short time t / and (b) final 
frequency cj/, where the parameters are the same as Fig. [2] 
polynomial trajectory (dotted blue line), A A relation (solid 
red line) . 



and the relation 



AHo tf > H arccos 



V^oUJf} 1 / 4 

(wo + Uf) 



1/2 



(23) 



In the Uf — > limit one finds again the time-energy 
uncertainty relation for two orthogonal states, that is, 
AHo > h/(4tf), independent of uif. This bound, al- 
though certainly correct, is not tight and does not cap- 
ture the actual dependences found for time averaged 
standard deviations, which in fact scale on utf and tf 
in the same way as the corresponding time averaged en- 
ergies of the previous section. Similar to Fig. [5J Fig. [S] 
makes the exponents explicit. We have used the stan- 
dard deviation AH n (t) = [(*„|iJ 2 (t)|*„) - E*] 1 / 2 for 
the n-th expanding mode, which takes the form 



AH n (t) 



V2(r 



n ■ 



4w 

,2 



L0\t)b' 



b 2 



-1 1/2 



b 2 



(24) 
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and its time average 



AH n = - / AH n {t)dt. 
tf Jo 



(25) 



The dependence of AHq ontf and w/ as they approach 
zero independently is summarized by the scaling expo- 
nents. In the limit of tf — !• 0, AH a cx A. Figure [S] 
shows that 6 = — 2 for the calculated standard deviation, 
whereas AA provides S = — 1. Similarly, as ujf — > 0, 
AHq oc uij. We find 6 = — 1 in the calculated standard 
deviation versus 5 = from the AA relation. 



V. ADDING TERMS TO THE HAMILTONIAN 

Motivated by recent experimental realizations [i"3| . we 
have considered up to now simple processes in which 
the only external manipulation consists in shaping u)(t). 
Other possibilities exist in which the Hamiltonian is com- 
plemented with additional terms [!, Q . We shall analyze 
here the energy excitations for a Hamiltonian that re- 
sults from the transitionless inverse engineering algory- 
thm proposed by Berry |3j (The application to the time 
dependent oscillator was worked out in p.) We assume 
here that uj(t) remains positive, 



H = H + Hi, 

H = huj(t)(a\a t + 1/2), 



(26) 
(27) 

(28) 



H would drive the system without transitions along the 
states of the instantaneous basis of the time-dependent 
harmonic oscillator H. In particular for the n-th state, 



■Jo e„(.t')dt'\ 



(29) 



is an exact solution of the time-dependent Schrddinger 
equation with H . 

The subscript t in the Schrodinger-picture creation 
and annhilation operators above denotes their funda- 
mental time dependence, because of the changing fre- 
quency and eigenstates, not to be confused with the 
time-dependence of Heisenberg picture operators. Hi is 
related to the squeezing operator @, and may also be 
written as Hi = —j^(qf> + pq), so that H is still a gen- 
eralized harmonic oscillator quadratic in positions and 
momenta. The expectation values of H and H 2 for the 
n-th state are easily calculated, 



(0n\H\4>n) 



\H Z 



(30) 



> 2 +n + l) + e 2 n , (31) 



8w 2 

AH n ^ {{<t> n \H 2 \K) - 4) 1/2 - T^[2(n 2 + n + l)] 1 / 2 , 



from which we deduce 



and the time averages 

_ ?i(n + l/2) [*f 



u(t)dt, (33) 



AH n = — [2( 



4f 



/ 



+ l)]i/2/ W dt . (34) 
Jo u 



We shall first show with some specific examples that 
it is easy to find scalings which are in principle more 
favourable for implementing a fast transitionless process 
than the ones in the previous sections. For example, for a 
linear frequency ramp expansion u = ujq + (w/ — Wo)t/tf, 



= H(n + 1/2) 
which is independent of t / , and 



UJ + UJf 



AH n 



to 



-[2(n 2 



f 



l)] 1 / 2 ^ 



wo 

Uf 



(35) 



(36) 



(32) 



(in fact a general result for ill < 0) so, for a fixed AH n , 
tf grows only logarithmically as ujf — > 0. 

Note that some of the difficulties with high transient 
energies in the approaches which only control the time de- 
pendent frequency, due to the particle exploration of re- 
gions far away from the trap center, disappear here since 
the system evolves at all times along the instantaneous 
eigenstates without transitions. Clearly (<f> n \H\(f> n ) = e n 
and {(j) n \H 2 \4> n ) — ((t> n \H\(j) n ) 2 = 0, so that the stan- 
dard deviation (|36[) is entirely due to the complementary 
Hamiltonian Hi. 

The main, and so far important difficulty with this 
approach is that it is not clear how to implement H in 
practice [||. Hi involves a non-local interaction and the 
attempts to provide a quantum-optical realization have 
not yet succeeded. 



VI. DISCUSSION AND CONCLUSION 

We have studied the transient energy excitation in time 
dependent quantum harmonic oscillators engineered so 
that the level populations at a final time are the same 
as the initial populations. We have considered first sim- 
ple processes in which the only external manipulation 
consists in shaping uj(t). The populations of the instan- 
taneous levels at intermediate times are, however, not 
preserved, so the transient excitation should be under- 
stood and possibly controlled. We have obtained bounds, 
shown examples, and determined the dominant depen- 
dences, which are different from the ones in the Anandan- 
Aharonov relation (22[. 

In a realistic application the oscillator will not be per- 
fectly harmonic and it is natural to set some maximum 
value to the allowed excitation. Then the minimal time 
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required for a fast expansion scales with the final fre- 

— 1/2 

quency as tf oc uj ^ .As the velocity determining step 
in quantum refrigerator Otto cycles this implies a depen- 

3/2 

dence R oc T c of the cooling rate, which had been pre- 
viously conjectured to be a universal dependence char- 
acterizing the unattainability principle for any cooling 
cycle [20j . The present results provide strong support for 
the validity of this conjecture within the set of processes 
defined exclusively by time-dependent frequencies (with- 
out added terms in the Hamiltonian) , and call for further 
testing and study. 

In Sect. [V]we have seen that, at least at a formal level, 
one could design even faster processes by adding terms to 
the harmonic oscillator Hamiltonian, but their physical 



implementation remains challenging. 
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